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Abstract 

The ground state entropy of the 2D Ising spin glass with +1 and -1 bonds is studied for L x M 
square lattices with L < M and p = 0.5, where p is the fraction of negative bonds, using periodic 
and/or antiperiodic boundary conditions. From this we obtain the domain wall entropy as a 
function of L and M. It is found that for domain walls which run in the short, L direction, there 
are finite-size scaling functions which depend on the ratio M/L ds , where d$ = 1.22±0.01. When M 
is larger than L, very different scaling forms are found for odd L and even L. For the zero-energy 
domain walls, which occur when L is even, the probability distribution of domain wall entropy 
becomes highly singular, and apparently multifractal, as M/L ds becomes large. 

PACS numbers: 75.10.Nr, 75.40.Mg, 75.60.Ch, 05.50.+q 
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I. INTRODUCTION 



Recently, a series of unexpected results and conjectures^^ has provided a new per- 
spective on the nature of the two-dimensional (2D) Ising spin glass. These results indicate 
that for the case in which the bonds are chosen randomly to have values of ±J conformal 
invariance^ and gauge invariance^ combine to create a set of remarkable properties. In this 
work we will present more such remarkable properties of the model, and find a numerical 
result which seems to rule out some of these proposals. 

The Hamiltonian of the Edwards- Anderson spin-glass model^ for Ising spins is 

H = - JijViVj , (1) 

<y> 

where each spin Oi is a dynamical variable which has two allowed states, +1 and — 1. The 
(ij) indicates a sum over nearest neighbors on a simple square lattice of size L x M. We 
choose each bond Jy to be an independent identically distributed quenched random variable, 
with the probability distribution 

P(J^ = P 5(J tJ + 1) + (1 - p)8(J ij - 1) , (2) 

so that we actually set J = 1, as usual. Thus p is these concentration of antiferromagnetic 
bonds, and (1 — p) is the concentration of ferromagnetic bonds. In this work we will study 
the case p = 0.5, for which the average of P{Jij) is zero. 

In two dimensions (2D), the spin-glass phase is not stable at finite temperature. Because 
of this, it is necessary to treat cases with continuous distributions of energies (CDE) and 
cases with quantized distributions of energies (QDE) separately.-^ 

Amoruso, Hartmann, Hastings and Moore^ have proposed that in 2D there is a relation 

where ds is the fractal dimension of domain walls, and 9e is the exponent which characterizes 
the scaling of the domain wall energy, E^w, with size. For the CDE case, the existing 
numerical estimates^- of ds and 6 E satisfy Eqn. (3). 

For the QDE case, it is known that 6% = 0. 9,11 Using Eqn. (3) then gives ds = 1.25. 
The derivation of Eqn. (3) assumes that the critical exponent r] for the scaling of spin-glass 
correlations is equal to zero, however. This appears to fail in the QDE case.— 



As pointed out by Wang, Harrington and Preskill, 15 domain walls of zero energy which 
cross the entire sample play a special role when the boundary conditions are periodic and/or 
antiperiodic in both directions and the energy is quantized. In the work presented here, we 
will find that the properties of these E^ w = domain walls are very special indeed. 

II. NUMERICAL RESULTS FOR GROUND STATE ENTROPY 

We will analyze data for the domain wall entropy, Sd w , for the ground states (GS) of 
2D Ising spin glasses obtained using a slightly modified version of the computer program of 
Gallucio, Loebl and Vondrak,— ^ which is based on the Pfaffian method. The Pfaffians are 
calculated using a fast exact integer arithmetic procedure, coded in C++. Thus, there is 
no roundoff error in the calculation until the double precision logarithm is taken to obtain 
Sdw- This extended precision is essential, in order to obtain meaningful results for entropy 
differences for large values of L and M. An earlier version of this calculation^ was limited 
to L x L lattices. 

We define the GS entropy to be the natural logarithm of the number of ground states. 
For each sample the GS energy and GS entropy were calculated for the four combinations of 
periodic (P) and antiperiodic (A) toroidal boundary conditions along each of the two axes 
of the square lattice. We will refer to these as PP, PA, AP and AA. The computer program 
treats the two lattice directions on an equal footing. Therefore we assume, without loss of 
generality, that L < M. 

For each lattice size L x M which was studied, 500 samples of the random bonds were 
used to calculate statistical averages for quantities of interest. In Fig. 1 we show the average 
GS entropy [Sq(L,M)] per spin, where the brackets [ ] indicate an average over random 
samples of the J^-, for a large number of sizes L x M. The values of M shown here are chosen 
to be 2n + 1, where n is 0, 1, 2 ... . Thus for these lattices L and M are either both odd or 
both even. We see that the behaviors for odd L and even L are distinct. When the aspect 
ratio, R = M/L becomes large, [S (L, M)\ for the lattices with odd L approaches S {oo, oo) 
from below, while for even L it approaches this limit from above. R functions as a control 
parameter which takes us from the L x L systems, for which [So] falls on a single curve for 
both odd L and even L,— to the L x oo systems, for which [So] has distinct behaviors for 
odd L and even L. 
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FIG. 1: (color online) Average GS entropy per spin, [Sq]/LM vs. 1/LM, for LxM lattices. The 
error bars indicate one standard deviation statistical errors. 

We define domain walls for the spin glass as it was done in the seminal work of McMillan.— 
We look at differences between two samples with the same set of bonds, and the same 
boundary conditions in one direction, but different boundary conditions in the other direc- 
tion. Thus, for each set of bonds we obtain domain wall data from the four pairs (PP,PA), 
(PP,AP), (AA,PA) and (AA,AP). For each size L x M we have 1000 data points for the 
short (horizontal) direction, and another 1000 data points for the long (vertical) direction. 

The domain-wall renormalization group 19 is based on the idea that we are studying an 
effective coupling constant which is changing with L and M. For the CDE case- we can 
use the domain wall energy, E^w, which is defined to be the change in the GS energy when 
the boundary condition is changed along one direction from P to A (or vice versa), with the 
boundary condition in the other direction remaining fixed, as the coupling constant. For 
the QDE case, what we need to do is a slight generalization of this idea. We should think 
of the coupling constant as the free energy at some infinitesimal temperature. When we do 
this, the entropy contributes to the coupling constant. 

The domain wall entropy, Sd w , is defined^ to be the change in S when the boundary 
condition is changed along one direction from P to A (or vice versa), with the boundary 
condition in the other direction remaining fixed. As long as E^ w > 0, the two boundary 
conditions which we are comparing are not on an equal footing. At a fixed aspect ratio, 
[Sdw] is expected to increase as a positive power of L for any E dw > 0. Therefore, these 
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FIG. 2: (color online) Finite-size scaling function for [\Sd w (L, M)\] vs. M/L 1,22 for Ed w = domain 
walls which run in the L direction. The y-axis is scaled logarithmically. 

coupling constants must eventually, at large enough L, be controlled by [Sdw] for any T > 0. 
Of course, the value of L which is needed for this to happen depends in T. 

As Wang, Harrington and Preskill 15 express the situation, an Ed w > domain wall 
does not destroy the topological long-range order. However, in the Edw = case the two 
boundary conditions are on an equal footing, and the topological order is destroyed. The 
probability distribution of Sd w for the cases where Ed w = should be symmetric about 0, 
and our statistics are consistent with this. Therefore the Edw — class of domain walls can 
be expected to behave in a special way. 

It is important to realize that the meaning of a domain wall is very different when Sq 
is positive, as in the model we study here, as compared to the typical case of a doubly 
degenerate ground state. In the case of two-fold degeneracy one can identify a line of bonds 
which forms a boundary between regions of spins belonging to the two different ground 
states. It is not possible, in general, to do that when there are many ground states. Despite 
this, we continue to use the term "domain wall". 

III. S dw FOR EVEN L 

When L is even, the energy difference, Edw, for any pair of GS for which the boundary 
conditions are changed in the horizontal direction, with the boundary conditions in the ver- 
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tical direction remaining fixed, must be a multiple of 4. When L is odd and the boundary 
conditions are changed in this way, Ed w is 4n + 2, where n is an integer.-^ Equivalent state- 
ments are true for odd and even M, with the roles of the horizontal and vertical boundary 
conditions interchanged. The sign of E^, w for a McMillan pair is essentially arbitrary for 
p — 1/2. Thus we can, without loss of generality, choose all of the domain- wall energies to 
be non-negative. 

The probability distribution for Ed w is a strong function of the aspect ratio of the lattice. 
When L is even the probability that E dw ^ goes exponentially to 0, as a function of R.— 
Similarly when L is odd we find that the probability that Ed w ^ 2 goes exponentially to 
zero as R increases. However, the difference between even L and odd L for large R turns 
out to be profound. 

For the E dw = case it is convenient to study \Sd w \, since the distribution is symmetric 
about zero. In Fig. 2 we display a scaling function for the behavior of [IS^I] versus M/L ds , 
for domain walls which run in the (L) direction, with 

d s = 1.22 ±0.01, (4) 

where the error estimate is a one-standard-deviation statistical error. This estimate of ds is 
obtained from comparing the behavior of the L = 8 data and the L = 12 data as a function 
of M, and requiring that the slopes in Fig. 2 should be identical. Thus we find that the true 
value of ds is probably less than 1.25, and the possibility that ds is the same for the QDE 
case as the CDE case,- for which ds is approximately 1.27,— appears to be ruled out. 

The naive choice of scaling variable, R = M/L, fails to produce a satisfactory data 
collapse for data with differing values of L. L ds is the effective length of a domain wall for 
this model. It is intuitively reasonable that the domain wall entropy should be proportional 
to the length of the domain wall, rather than L. Although such data are not displayed here, 
we have found that when M > L the results for odd M with even L also fall on the scaling 
curve shown in Fig. 2. 

Recently, Melchert and Hartmann 2 ^ have attempted to calculate the scaling exponent for 
the length of domain walls for this model directly. They were unable to determine a precise 
value, however, due to the high ground state degeneracy and the fact that their algorithm 
does not select a ground state randomly. Similar results for hexagonal lattices have been 
given by Weigel and Johnston. 24 
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FIG. 3: Histogram of the probability distribution of [15"^ (6, 66) |] for E& w = domain walls which 
run in the L direction. The y-axis is scaled logarithmically. 

It is remarkable that [IS^I] is falling exponentially as a function of the variable M/L ds . 
There is no L-dependent scale factor for the y-axis. This is the total \Sd w \ for the entire 
L x M lattice. This behavior indicates that the zero-energy domain walls which encircle the 
lattice in the short direction must become strongly correlated as R becomes large. Because 
the Hamiltonian does not contain any explicit long-range interactions, the mechanism by 
which this occurs is not trivial to understand. 

As M/L ds increases, the distribution of [ | -S 1 ^ | ] becomes increasingly singular. The last 
point shown in Fig. 2 is above the trend because it is dominated by a single data point. To 
make this issue concrete, in Fig. 3 we show the histogram of the probability distribution 
for this point, 15^(6, 66)|. The domination of the mean of this distribution by the point at 
the far right of the histogram is obvious by inspection. To demonstrate that this is typical 
behavior, we show the corresponding histogram for \Sd w (12, 132) |, which is not off the trend 
line, in Fig. 4. 

As the reader may have already noticed, no statistical error bars are given for the data 
points in Fig. 2. In order to give a meaningful estimate of such statistical errors for these 
probability distributions, we would need to know what the analytical forms of the probability 
distributions are. We do not have this information. The statistical error estimate for ds 
comes from fitting the observed fluctuations of the data points from the trend line. For this 
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FIG. 4: Histogram of the probability distribution of [15^(12,132)1] for Ed w = domain walls 
which run in the L direction. The y-axis is scaled logarithmically. 
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FIG. 5: (color online) Finite-size scaling function for [\f\Sd w {L, M)\] 2 vs. M/L 1 - 22 for Ed w = 
domain walls which run in the L direction. The y-axis is scaled logarithmically. 

we do not need to know the statistical errors of the individual data points. 

Because the probability distributions of \Sd w \ become so singular in the limit of large 
R, they appear to be multifractal^. One way of seeing this is to calculate the fractional 
moments 

\S dw {L,M)\ q = [\S dw (L,M)\^]\ (5) 



S 



for q= 1, 2, 3 ... when R > 1 for the Edw = domain walls which run in the L direc- 
tion. If \Sdw\ was controlled by a single length scale, the slopes of the scaling functions 
for \Sd w (L, M)\ q , analogous to the q = 1 case shown in Fig. 2, would be identical (within 
statistical errors). When one does this calculation, however, one finds that the rate at which 
\Sdw(L, M)\ q decays exponentially to zero as M/L ds increases is an increasing function of q. 
The results for q = 2 are shown in Fig. 5. It is clear that the data points are falling faster 
as M/L 1,22 increases for q = 2 than they do for q — 1. (Note that the scales on the axes of 
Fig. 5 are different from those of Fig. 2.) This trend continues for larger values of q. 

The probability distributions for \Sd w \ of Edw — domain walls are thus inconsistent with 
any single-length scaling hypothesis. The author believes that this behavior arises from the 
fact that the eigenvectors of the susceptibility matrix must be Anderson-localized in 2D. 26 
Mult if r act al behavior is a natural result of disorder-induced localization. 25 

In an earlier work,- the author calculated the finite-size scaling behavior of \Sd w \ for 
Edw — domain walls on L x L square lattices. It was found there that the exponent 9 s for 
the scaling of the width of the probability distributions for \Sd w \ depends on Ed w = 0. For 
Edw = domain walls 

e s (E d w = 0) = 0.500 ± 0.020 . (6) 

According to the droplet model, one expects ds = 26s-— The value of ds which we have 
found here is clearly inconsistent with this relation for Edw — 0. The reason for this is the 
special symmetry of the Edw = domain walls, as discussed in the earlier work.- 

IV. S dw FOR ODD L 

For odd L with the boundary conditions we are using, Ed w cannot be zero for domain 
walls which run in the L direction. When Edw is n °t zero, the relative signs of Edw and Sdw 
are not arbitrary. Having chosen Edw to be positive, we then find that Sdw is also usually 
positive.- In Fig. 6 we show the finite-size scaling of [Sd w {L, M)\ for the Ed w = 2 domain 
walls which run in the L direction. In order to make the data for different values of L fall on 
a common curve, it is necessary to include a y-&xis scale factor of L~ l l 2 . The scaling variable 
again appears to be M/L ds , with the same value of ds as before, although our uncertainty 
in the value of ds is larger for odd L. The behavior seems to be a power-law increase, with 
an slope close to 0.30. The contrast of this behavior with the exponential decrease which 
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FIG. 6: (color online) Finite-size scaling function for [Sd w (L, M)]L vs. M/L for E t 
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domain walls which run in the L direction, log-log plot. 
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FIG. 7: (color online) Finite-size scaling function for the standard deviation, sd(Sd w {L-,M))L 0,67 
vs. M/L 1 ' 22 for Ed w = 2 domain walls which run in the L direction, log-log plot. 

we found for the Ea w = walls could hardly be greater. 

In Fig. 7 we show the finite-size scaling behavior for the widths of the same Ed w = 2 
Sd w (L,M) distributions, as parameterized by their standard deviations. The y-axis scale 
factor for data collapse appears to be L~ 2//3 in this case. The behavior again appears to be 
a power-law in M/L ds . The slope in this case is approximately —0.30. Thus the widths 
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of these Sd w {L,M) distributions should eventually become small compared to their average 
values as R increases. 

For domain walls that run in the long, M direction, the probability that Ed w = goes to 
zero as R — > oo.— We do not expect any anomalous behavior for the "long" domain walls. 
All of the results found here are consistent with the conclusion 3 that in this model there 
seem to be two distinct classes of domain walls, the Ed w = domain walls and the Edw > 
domain walls. 

The alert reader has noticed that the sizes of the lattices for which we have data when 
L is odd are more limited than in the even L case. This is because the Vondrak code runs 
about a factor of four faster if both L and M are even, due to symmetry properties of the 
Pfaffians. 

V. SUMMARY 

We have studied the statistics of domain walls for ground states of the 2D Ising spin 
glass with +1 and —1 bonds for L x M square lattices with p = 0.5, where p is the fraction 
of negative bonds, using periodic and/or antiperiodic boundary conditions, for both even 
and odd L and M, where L < M. The probability distributions of domain wall entropy, 
Sdw(L, M), are found to depend strongly on E dw , and therefore on whether L is odd or 
even. Finite-size scaling forms are found which are functions of the variable M/L ds , where 
d s = 1.22 ± 0.01. When the aspect ratio becomes large, the distribution of Sd w for zero- 
energy domain walls which encircle the lattice in the short direction becomes multifractal. 
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